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Abstract 

o . 

^ ^ Most of the special functions of mathematical physics are connected with the representation 

^\ ■ of Lie groups. The action of elements D of the associated Lie algebras as linear differential 

operators gives relations among the functions in a class, for example, their differential recurrence 
\l , relations. In this paper, we define fractional generalizations of these operators in the context 

of Lie theory, determine their formal properties, and illustrate their use in obtaining interesting 
relations among the functions. We restrict our attention here to the Euclidean group E(2) and 
the Bessel functions. We show that the two-variable fractional operator relations lead directly 
to integral representations for the Bessel functions, reproduce known fractional integrals for 
f~>) I those functions when reduced to one variable, and contribute to a coherent understanding of 

. the connection of many properties of the functions to the underlying group structure. We extend 

the analysis to the associated Legendre functions in a following paper. 

(N : 
o 

I ■ Most of the classical special functions are connected vi^ith the representation of Lie groups |[ ^, |^ , 

I and appear as factors in multivariable functions on which the action of an associated Lie algebra is 

■ realized by linear differential operators. Many of the properties of the special functions are easily 
^ I understood in this context. For example, the differential equations for the special functions are 

■ connected with the Casimir operator of the associated groups. The actions of appropriate elements 
. !^ I 13 of the Lie algebra lead, when reduced to a single variable, to the standard differential recurrence 

■ relations for the functions, while the action of group elements e~*^ can be interpreted in terms 
H I of generalized generating functions when expressed using a Taylor series expansion in the group 

parameter t. Numerous examples are given in In the present paper, we will define fractional 

generalizations Z)^ of the D's in the context of Lie theory, determine their formal properties, and 
illustrate their usefulness in obtaining further interesting relations among the functions, including 
integral representations for the functions. Most of the specific results have been derived historically 
in other ways, but are unified here in a group setting. 

Two example of fractional operators in a single variable are provided by the fractional integrals 
of Riemann and Weyl Chap. 13. These give a useful way of changing the indices (degree or order) 
of the classical orthogonal functions (Jacobi, Gegenbauer, Legendre, Laguerre, Bessel, and Hermite 
functions). An example is Sonine's first integral for the Bessel functions 12.11.(1), 

x'+f^J^+^ix) = I ' •^-(^Sin0)cos2^-i0sin^+i0d0 (1) 



1 Introduction 
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f+\j4t){x^ ^t^)''-Ut. (2) 



The expression on the second Une is equivalent to the Riemann fractional integral Rfj^x'^^'^ J^{2^/x), 
where the integral operator i?^ is defined in general by Chap. 13, as 

{Rc.f){x) r f{t)(x - tr-' dt. (3) 



r(a) Jo 

Thus, with the replacement of t by 2\/t and x by 2a;, (||) becomes 

^(.+M)/2j^^^(2V^) = -i- r t-'^M2Vt){x-tY-Ut (4) 

= R^,x''/^J^{2^). 

A number of similar results are known for other special functions, for example, the relation 

2^-M(i _ x^)-'^^-^'^/^P^-^'{x) = i?^2^(l - x^)-^/^P^{x) (5) 

for the associated Legendre functions with ReA < 1, Re/^ > 0, ||], 13.1(54). Askey 0, Chap. 3, 
summarizes a number of results and gives some applications. 

Other results are known with respect to the Weyl fractional integral W^^ Chap. 13, defined 

by 

{Wo.I){x) = — - / f{t){t - xr-' dt. (6) 
r(a) Jx 

Thus, from §, 13.2(59), 

^-(.-M)/2j^^_^(2^) ^ w^x-''K,{2y/^) (7) 

where K^, is the hyperbolic Bessel or MacDonald function. 

The simplicity of the results noted, and of many similar results is striking. The effect of 
the fractional integration is simply to change the indices on the special functions, while retaining 
the original functional form. There does not appear to be a systematic approach to the derivation 
of these results in the literature. Their form suggests that they must be associated with fractional 
generalizations of the stepping operators in the associated Lie algebra. In particular, the differential 
recurrence relations for the special functions are schematically of the form DFa^... — ci^Q±i,..., 
where D is a linear differential operator and the indices a label the functions in a realization of 
the algebra. This suggests that shifts of the indices by arbitrary amounts could be effected using 
fractional operators defined in analogy to the single- variable fractional derivatives defined in j|] . 
This is the case, as we will see. The fractional integrals given above are related, and simply give the 
action of the inverse multivariable operators D~'^ when reduced to a single variable. 

We will define the fractional operators D^^ in the context of Lie theory and explore their general 
properties in the following section. We will then apply the results in a number of group settings in this 
and following papers to obtain generalized fractional- integral- type relations of the form i^Q+^,... = 
ND^^Fa^... for the special functions. Some are apparently new. We find that, with appropriate choices 
for the input functions, the fractional relations lead directly to known integral representations for 
the special functions, providing a group-theoretical setting for the latter. 

In the present paper, we will restrict our attention to the development of our methods, and to 
applications to the Bessel functions. Our treatment is not exhaustive in either the theory or the 
applications considered. 
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2 Fractional operators 



We will suppose that we have a Lie algebra which corresponds to one of the classical Lie groups, 
and is realized by the action of a set of linear differential operators {D{w,dnj)} in a collection of 
variables w acting on an appropriate set of functions. The exponentials e~*^ defined by Taylor series 
expansion in the group parameter t are elements of the Lie group taken to act on an appropriate 
class of functions F. We will assume that the group action e~*^F can be defined for all t, and will 
define a Weyl-type fractional operator as an integral over group elements by 

D^iw, d^)Fiw) = 2^e'-''r(/. + 1) / dt Fiw). (8) 

The contour Cw = (oo, 0+, oo) in the complex t plane runs in from infinity, circles t — Q in the 
positive sense, and runs back to infinity. To define phases, we take the integrand as cut along the 
positive real axis with the phase of t taken as zero on the upper edge of the cut. The direction of 
the contour at infinity must be such that the integral converges. The expression above would be an 
identity for D a positive constant. Here, however, D{w,dw) is an operator which acts on functions 
F of the collection of variables w, and the existence of the integral depends on the functions as well 
as the contour. 

Alternatively, D^y can be defined as 



D^lF = — ^ D" / dt — F, (9) 

where Re fj, < n and endpoint terms are assumed to vanish in the partial integrations which connect 
the two expressions. 

It is straightforward using this expression to show that the fractional operators have the expected 
algebraic properties. Thus, for Re/i < 0, Kcv < 0, and Rc(// + < 0, 

1 j-OO j-OO 1 1 
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-M)r(-j^) 
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dv^—r- / dt't'-^-\l~t')-"~' (10) 

^'^+"^+' Jo ^ ^ 



CO p—vD 

dv^^=D>^''F 



r(-/l - Jo 



Exponents therefore add as we would expect, and the fractional operators of different orders com- 
mute, 

D^D^w = D'^D^ = D^''^ [D^,D^w]=0- (H) 

The result extends throug h pj to general /i, v for which the fractional operators are defined. By 
converting the integral in [10| back to a contour integral before taking the limit v — > — /z, we find 
also that D'^D^ = 1 where 1 is the unit operator, so is the inverse of D^^ as implied by the 
group operations. 

The fractional operator D!^ can also be defined in terms of the action of a generalized Weyl 
fractional integral in the parameter x on the group operator g~^D{w,d^) ^ .^jjj define W-^, 
for general fi as 

1 ..^MT.... , ,^ /■ fit) 



W.M = -e-r(, + 1) J^^ dt dt, (12) 
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where Cx is the contour (oo,a;+,oo). This definition is equivalent to [s] for Re/x < 0. The action of 
W-^ on / is just that of a fractional derivative, 

{W^^f){x) = {-d/dxYfix), (13) 

a result which is obvious for /i an integer so that the integration contour can be closed. In general, 
W-^ gives the inverse of thought of as a repeated integral, W-^W^ = 1. 

Dy^ can now be defined formally through the action of the fractional derivative {—d/dx^ on 
e"^^, {-d/dxYe-^'= = D'^e"^^. Multiplication by e^^ then gives = {W-^,e'^'-'° F). 

This relation is easily checked by using J{t) — e~^^ F in |lj and changing the integration variable 
from t to t — X. We find that 

1 P p"*^ 

D'^F = e^^{W^^e-^'')F = —e^^m^l + 1) dt —^F (14) 

in agreement with |8[ That is, D^^F — e^^ W-^{e~'-^^ F ) where W-f^ acts on the group pa- 
rameter X and D{w,dw) acts on F{w). The integrals in ^ and ^ can also be identified di- 
rectly as {—d/dx)'^e~^^\x=Q = D'^. The inverse of the fractional operator D'^F is F>^F = 
e^°Wf,{e-''^F). 

We can define a second Riemann-type fractional operator by replacing the Weyl fractional integral 
by a Riemann fractional integral and noting the correspondence of R-a to (d/dx)", R-af{x) = 
(d/dxYfix). Thus, taking taking f{t) = e*^F and a = n — Re/i > in ^ and following the 
construction above, we find 

D^F = e-^-°L>"(i?„_^e^^)F, 0<n- Refi. (15) 

By changing the integration variable from t to a; — t in ^ we then obtain the analog of 

1 i-x{w) tD 

{D'tF){x) = — -D" dt —F, (16) 

where we have noted the dependence of the final result on the value x{'w) of the group parameter t 
at the endpoint of the integration. As indicated, this will depend on the values of the variables w 
in F. 

By going to a contour integral to handle the possible singularity at the lower limit of integration, 
we can write D^F in the more general form 



^5i^-ie-r(, + i)/^^d.|^F, (17) 



where Cji is the contour Cr — {x(w), 0-I-, x{w)). 

As we will see explicitly in later applications, the endpoint x{w) of the contour must be chosen 
such that D^F satisfies a differential equation determined by the Casimir operators of the Lie 
algebra. This will require that a differential expression related to e~*^F vanishes for t — x{w) for 
the given values of the variables w in F (see, for example, ^5| ). 

The product of two Riemann fractional operators is given in the simple case Re /i. Re < by 



dt du^^ —e(x-t~u). (18) 



D'^{D'^F ) will satisfy the expected differential equation for D^G provided t + u = x on the boundary 
of the region of integration, a condition is enforced in Eqby the unit step function 9{z), 6{z) — 1 for 
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z > and 9{z) = for z < 0. For an explicit example, see §3.4.1. The integral can be evaluated by 
shifting to ti = t + u as a new integration variable and identifying the remaining integral with a beta 
function as in |lO|. The result is equal to D^'^ F. We therefore obtain the multiplication relation 
^R^R = ^R^R ~ -^fl^" derived earlier for the Weyl fractional operators. This can be generalized 
to the contour integral representation 

Which expression for D'^ is appropriate in a particular setting, Weyl or Riemann, will depend on 
D and F. We will therefore simply denote the fractional operator as for formal purposes, and 
only specify the expression to be used in in particular applications. The key restrictions will be the 
existence of a finite value of the group parameter t = x{w) such that e"^^ F = in the Riemann 
case, and convergence of the integral for i ^ oo in the Weyl case. 



3 Bessel functions and E(2) 
3.1 Algebraic considerations 

As a first application of the fractional operators, we will consider the Bessel functions which we 
will denote generically as Bessel functions appear naturally in representations of E(2), the 

Euclidean group in two dimensions, and of E(l,l), the Poincare group in two dimensions 1^, Both 
groups are real forms of S0(2,C), and the two are related to eachother through the Weyl unitarity 
trick Since we are not concerned with unitary representations of the groups, it will be sufficient 
for our purposes to consider only the algebra of E(2). 

The Lie algebra of E(2) is generated by three operators Pi, P2, J3 with the Lie products or 
commutation relations 

[^1,^^21=0, [J3,Pl]=P2, [J3,P2] = -Pl- (19) 

There is one invariant operator, namely + P2, which commutes with all the generators. 

The algebra can be realized by the action of differential operators on functions / of the coordinates 
{xi, X2) in the Euclidean plane. Pi and P2 correspond to the translation operators 

Pi=du P2 = d2 (20) 

and J3, to a rotation in the plane, 

J3 = -xid2 + X2di. (21) 

The condition that the invariant operator P^ + P| be constant on the functions / gives the 
Helmholtz equation {P^ + = —k^f- In polar coordinates x, (f) this becomes the differential 

equation 

m^ + ^l + ^V^[|^4^ + ;^^ + '.-]/^0. (22) 

We can take fc^ = 1 by a scaling of the coordinates, and will do so. The rotation operator J3 ~ —d^ 
commutes with the Helmholtz operator and may also be taken to have a constant value —iiy on the 
functions. The functions / in this realization of E(2) are then of the form 

Mx,<j>)=e'''^Z,{x), (p2 + p| + i)/ = o, hfu = -iyU (23) 

and involve Bessel functions Z^, of order v. 

It is useful to change from the antihcrmitian operator J3 to the hermitian operator iJ^, and to 
introduce operators 

P+ = -Pi - iP2. P-=Pi- iP2 (24) 

with the commutation relations 

[P+,P-] = 0, [iJ3,P±] = ±P±. (25) 
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The last relations imply that if fi, is a solution of the Helmholtz equation with the index v, then 
P±.fv is a solution with index v ±1, 

iMP±U) = P±{^J3 ± l)U = (i^ ± l)iP±U)- (26) 

P± therefore act as stepping operators on the index. 
The operators are given explicitly by 



where t = e*"*. In terms of that variable, i J3 — tdt- The Helmholtz operator is simply P+P- + 1. 

From ^ the action of P± on the functions f^{x,t) — Zi,{x) must give constant multiples of 
t^^^ Zy±i{x). The constants of proportionality for the different Bessel functions are easily determined 
to be unity by using the behavior of the functions for x 0,oo. We therefore have the stepping 
relations 

P±t''Z,{x)=f'^'Z,±,ix), (29) 

which reduce to 

( T ^ + ^) Z,{x) ^ Z„±i{x) (30) 

once the t dependence is factored out. The latter are just the differential recurrence relations for 
the Bessel functions 7.2.8. 

The relations in 29 suggest that 

F£rZ,(x) =r±^Z,±^(x) (31) 
for appropriately defined fractional operators such as the Weyl operators 

P± = 7^e"''r(Ai + 1) / du (32) 

It is easily established that these operators have the expected properties. First, P!^ commute with 
the Helmholtz operator P+P- + 1, so transform solutions of the Helmholtz equation to solutions. 
Further, from the relation 

[iJ3,P£] = ±nP£, (33) 

we find that 

[zJ3,e-"^±] = ±^^L„p£ (34) 

n=0 

= ±u^e-^P±, (35) 
du 

hence, after a partial integration in that 

[zJ3,P£] =±mP£. (36) 
The action of P^ on a solution therefore gives another solution with the index v changed to vifJ,, 

lUP^U) = P£(*J3 ± m)/. = (^ ± ^^){P^f)■ (37) 
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This relation does note show directly that P±fu — fu±in but only that P±fi, is at most a linear 
combination of the two independent solutions of the Helmholtz equation with orders v ± ^l. If the 
independent solutions are taken as the Hankel functions, the observation that the operators P'^ do 
not change the distinct asymptotic behaviours of those functions for |x| — > oo shows, in fact, that 
Pl^-fu — N{h', m)/i^±a'- The constant of proportionality will be found later by direct calculation to be 
unity, as in establishing the validity of |3l|. 

We can also define a fractional operator (iJ^)^, and find after a brief calculation using the analog 
of |3^ that 

{^J3)^U = y^U- (38) 

[iJs)^ again satisfies the multiplication rule, [iJ^j^iiJsY = {iJi)^^^- 
The formal algebraic structure is completed by the relations 

(zJa) V£ = P£(*J3 ± ^Ji)\ P^iiJ^)^ = (zJs T (39) 

These can be derived using the Baker-Hausdorff expansion of e^Be^^ as a series of n-fold commu- 
tators, 

oo _j 

e^Be-^ = B + ^ - [^, [^, . . . [A, B] . . .]]. (40) 

n=l 

Thus, choosing as the exponential in the definition of {iJ:})^, = e"**"^^, B as P^^, and using ^ 
to evaluate the repeated commutators, we find that 



P± + y] ^—f- [*-^3, *J3, . . . [tJs, P^] ■ ■ ■]] e-'"-^ = p^e-*(*./3±p). (41) 



n=l 



The first of the relations ^ then follows upon integration using the analog of 32 , Application of this 
operator to a solution of the Helmholtz equation gives 

{iJzfP^U = PU^Js ± = ('^ ± f^)^P^U (42) 

The second of the relations ^ can be derived similarly. The complete algebraic structure defined 
by (i^a)^, P±, the multiplication rules, and |39| is infinite, and has not been investigated except as 
applied to solutions of the Helmholtz equation. 

3.2 Action of the group operators 

The action of the exponential operators 6""''^='= — g^'^Pi+^P^) jg easily determined and well known. 
Pi and P2 commute, and the exponentials e"^^ and e°'^^ induce translations of the coordinates xi, 
X2 with e'^^^xi = xi+a and e°'^^X2 — X2 + a. Thus, acting on functions analytic in the neighborhood 

of {xi,X2) 

e-^P+F{xi,X2) = e^P'e'^P'F{xi,X2) = F{xi + u,X2 + lu). (43) 
Applying this result to the functions Ju = e^^"^ Z^{x) written in rectangular coordinates, we find that 

\Xi-lX2j V 

= fx" {x^ + 2uxty''''^ (v/a;2 + 2uxt) , (44) 



where x = ^/x\+x^ and t — e*"^ = {xi + 1x2) /x. A similar calculation gives 



e--/. = ( - 1 [x^-^-:f]'\.\ \lx^ - ^ ) . (45) 



t\ / , 2ux 
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n 

n=0 



We can also calculate directly in polar coordinates, a method which will be useful later. Thus, 
noting that P+t"x'^ = {t/x){tdt — xdx)t'^x'^ = and using ^ we find that 

P;'rZ,(x) = (-2)"t''+"x'^+" (^^^ {x--'Z4x)) . (46) 

The formal Taylor series expansion of 6^"^^+ then gives 

= fx" {x^ + ^uxty"''^ Z^ (V^M^^^^ ^ (47) 

where we have identified the exponential in the penultimate line as a translation operator. The 
result agrees with A similar calculation for e~"^" reproduces |4^. 

Direct evaluations of e^'^^^ji, using the Taylor series for the exponentials and the relations 
P^t"Z^{x) = t''^^Zy±n{x), |2^, give the generating functions 

fx" {x^ + 2uxty'"''^ Z^ {^s/x-^ + 2uxt^ = i^^t''+"Z^+„(a;) (48) 

n=0 ^' 

and 

'^'^(^^ ( Jx^ - '-^] ^ ± t^t"^-Z.-Ax). (49) 



, , X . , ,, , y 

x J \ t I \V t I ^ — ' n 

^ ^ ^ \ / n— 

These equations give generalizations of Lommel's expansions for the Bessel functions Q|, §5.22. 
Thus, taking x = ^/z, ut = h/2x — h/2s/z in |4^, and choosing Zy as the ordinary Bessel function 
J^, we obtain §], 5.22(1), 



n=0 



(50) 

Similarly, for a; = and u/t = -h/2^, || gives @, 5.22(2), 

(z + hY'^MV^i) = y ^^z('^-")/V._„(yi). (51) 

^ — ' n! 

ri=0 

The remaining Lommel-type formulas in §5.22, follow from|4^ and |49| for different choices of Z^. 
The present development provides a group-theoretical derivation of these results. See also Weisner 
[ pO[ . Note the series ^ converges for |m| < |cc/2t|, and that in for \u\ < \xt/2\, that is, for 
sufficiently small values of the group parameter u. 

3.3 Weyl-type relations for Bessel functions 
3.3.1 Relations using 

The action of the Weyl-type operators Pj{^ on the Bessel functions is given by || and |4j, 
Ptt"Z,{x) - N{v,^i)t"+^^Zl^^{x) 

1 ,, ,, . ^ /" du 

~ 2TTi 



-t"x"e''"'T{fi+l) [ (x^ + 2uxt) Z^ (Vx^ + 2uxt \ . (52) 
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Cw is a contour {oo, 0+, oo) in the complex u plane with the direction of approach to oo to be taken 
such that the integral converges. This will depend on the function considered. 

Proceeding formally, we can extract the expected factor t^^^ from the integral by the change of 
variable v = 2uxt. We will also replace x by -y/z, with the result 



X 



'^\{y + xr'''^Z,{V^z). (53) 



This result can also be obtained directly from the differential recurrence relations ^ by replacing x 
by rewriting the resulting relation for P+ in the form 

- 2^ (z-^/^Z,{^,)) ^ (54) 

and determining the Weyl action of {—2d/dz)'^ on z~'^/'^ Zy{-s/z). 

The function z^^^"^ Z\{'s/z) satisfies the differential equation 9.1.53, 

Applying this operator with A = i/ + /i to the integral in converting the derivatives with respect 
to z to derivatives with respect to w, and using the differential equation for A = to eliminate the 
derivative- free term proportional to l/4z on the right hand side, we find that the result vanishes 
provided 

d { I d 



dv . 

(j^. dv y dv 



(i; + z)-'^/2^,(V^7TI)]| =0. (56) 



That is, the integral in |5^ gives a Bessel function or combination of functions with argument ^/z and 
order v + fi multiplied by provided the function in curly brackets vanishes at the endpoints 

of the integration. 

When Zi, is the Hankel function hI^\ the condition in ^ is satisfied for contours that run to oo 
in the upper half plane, avoiding the possible singularity at v — —z on the right. It also holds for a 
contour along the positive real axis for Re(/i + + f ) > 0. In either case, an asymptotic argument 

shows that the Bessel function Z^^^ given by the integral is in fact H^^^{^) with coefficient 
N{v, /^) = 1 as expected. Thus, 



X 



dv 



^^^-(« + z)-'^/^i?W(V^), (57) 



e > 0. Tracing the calculation back, we find that the original expression |5^ holds for for 
contours with < arg {xtu) < 27r as |w | ^ oo. 

A shift of the integration variable brings ^ to the form of a (generalized) Weyl fractional integral, 

^ J_2»e-(.r(p + 1) 

-'(oo,2 + ,oc) v*^ ^} 

The contour can be collapsed for Re/i < 0, and |5^ reduces to the known fractional integral ||], 
13.2(45). The latter can be written in the present notation as z-(''"^)/2i7iL'^(A/i) = P^^ z^" H^y\^/z), 
Re/x > 0. 
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(2) I 1 

Similar considerations for the choice = iJ^ ' in |52| show that that result holds for > 
arg{xtu) > 27r, and that 



z 



1 



2m 

-^(«+.r/^ij(^)(v^) (59) 

Cw 

where v runs to oo in the lower half plane, avoiding the possible singularity at w = ~z on the right. 
The result also holds for a contour along the positive real axis for Re(/i + + |) > 0. 

Combinations of and give the ordinary Bessel functions and the relations 

ZTTl 



% (t^ + z)""/' J.(V^), (60) 



Zm 

dv 

X 



- (z; + z)-''/'l^,(V^) (61) 



for Re(/^ + + |) > 0. The contours in these cases must be taken parallel to the real axis for 
V oo. The results reduce to the known fractional integrals 13.2(34) and 13.2(40) for Re/x < 0. 

If we increase the phase of z by tt and simultaneously rotate the contour Cw by tt in the positive 
sense in the expression p% the substitutions z = e"x, v = e^^u restore the original contour while 
replacing y/v + x by e^^^\/u + x. The definition of the MacDonald function K^, in terms of the 
Hankel function H^v^ ^ 

ir,(x) = ye"'^/2i7(i)(e»-/2a;), (62) 



then gives 



^-(.+;.)/2^ (yj) = -L2''e-^^r(/i + l) 

2-^1 



'^"-(u + x)-'^/2^,(V^7T^), (63) 



or, for Re/i < 0, 



in agreement with |^ or |5|, 13.2(59). 



3.3.2 Weyl-type relations from 

The action of the Weyl-type operators on the Bessel functions is given by || and |4£ 

1 ft 



2711 V X 
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Cw is again a contour (cxd, 0+, oo) in the complex u plane with the direction of approach to oo to be 
taken such that the integral converges. We will scale out the t dependence through the substitutions 
V = 2ux/t and x = y/z, and work with the reduced expression 

,i-^-f^)/2z (^^) ^ J_2>^e^^mp + 1) / -St(^ - vr/'Z,{^I^). (66) 

We will suppose initially that argz > 0. It is then possible for the choice Z,^ — H'iF' to rotate 
the integration contour into the lower half v plane. Then with v replaced by e^^'^v and z — v hy 
e'^'iv + z), 

.i^-.y^Hi^iV-z) = ^2-e^"-r(M + 1) J^^ -^{v + zY^'Hl'W^z), (67) 

where Cw is a contour (oo, 0+, oo) in the new variable v and — tt < arg z < n. By choosing Z^, = Hi 
and Imz < and rotating in the opposite sense, we obtain the second relation 

^(.-.)/2^(2) = ±2mp + 1) / -^{v + zy/'Hl'W^z), (68) 

also valid for — tt < argz < tt. These results can also be obtained by considering P^^V^ Z^- 
For Re/i < 0, the contours can be collapsed, and 

^(.-.)/2^(i^2)(^) = fni)^^"' [ + z)''/2i/(i'2)(V^), (69) 

where Imv ±oo for i?'^^ and H'^'^\ By considering the limiting behavior for Imw and 
combining the two functions, we obtain the relations 

- {v + zY^^J4^/v + z) 



= z^''-^)/^ [cos7r^J,_^(V^) + sin7r^i;„^(VJ)] , (70) 

2^ /-^ dv 



- {v + zy^^Y4V^) 



= z(''-A')/2 [cos7r^y,_^(V^) - sin7r^J,_^(V^)] . (71) 

These are equivalent to 13.2(35) and 13.2(39) and arc valid only for iRez^— | < Re/i < 0, with 
Im V —^ for Re w — > oo 



3.3.3 Weyl-type integral representations for Bessel functions 

We can use the results above to obtain integral representations for the Bessel functions. We begin 

[^^\{x) and _ffjy2 



with the observations that t^Z^{x) = ^ t^^"^ Zi/2{x), and that H^^}^{x) and h[^}^{x) are elemen- 



tary functions. 



The action of ^ can be reduced as above, and we will begin with the expression in |5^. This 
gives 

x-'^/'Hl!\^) = -L2M-ie^(.-i)^r(/. + i)|^ + 

= --L^e--r(M+i) / ^(. + x)-V2e^^/^. (73) 
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Replacing x hy x^, letting v ~ x'^ (t'^ — 1) , and removing a common factor of a; we obtain 

2Tn^\xJ 2 1+ (t2 _ 

This holds for general values of /i provided Ima;t oo for \t\ — > oo, and for — > +oo for Re/i > 
iO. The contour can be collapsed for Re/i < i giving the generalized Mehler-Sonine integral 
representation for H^^\x), ||] 6.13(1) 

The result satisfies the Bessel equation for Ixnxt oo for \t\ oo. 
A similar calculation for i?^ gives 



or, for Re/i < ^, 



(76) 



^ ^ ^ V^\x) r(i-/i) A (t2„i)M+i ^ ^ 

For x real and — ^ <Re/i< and |7^ can be combined to obtain the representations for 

and noted in 6.13(3) and (4), 

N 2 /2\^ 1 f"" , cosa;i 

A different set of integral representations can be obtained by considering the action of the inverse 
operator P;:^"* on t^^^ h[^/^\x), 

P^''''^t^/^H^^f^\x)=t-^H'^}f{x). (80) 

Using the relations 

(x) = e-^^i) (x), H^^lix) = ^42) (g^) 
and following the manipulations above, we obtain the integral representations 

ijW(.x) = '-^C'Ye-^^^^Til-t,) [ dt{t^~ir~h^^\ (82) 

TT V2/ 2 J(oo,l + ,oo) 

iif(.) = ^^f^Ynl-f,) f dtit'-ir-h-^^\ (83) 

TT V2/ 2 J(oo,l+,oc,) 

where, for convergence, t must approach 00 on the contours in ^ and ^ with Ima;i — > +00 and 
Im.Ti ~f —00, respectively. These expressions are equivalent to the representations 6.11(4) and 
6.11(5) in Q obtained from Hankel's representation for the Bessel functions.]^ Other results, for 
example, Schafli's integral for if^, |^] 6.15(4), can be obtained from these. See Watson ^. 

^Watson uses a different phase pnn ventinn in his 6.11(4) which is equivalent to replacing t — 1 in by e^'^*(t — 1). 
Watson's 6.11(5) is obtained from p3 by the substitutions t — 1 e"^(n + 1) and t + 1 ^ e^"^(u — 1). 
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3.4 Riemann-type relations for Bessel functions 
3.4.1 Relations for P£ 

For Riemann-type fractional operators, the roles of P'^ and are essentially reversed, and the 
relations apply to different Bessel functions. The action of the Riemann operator is given by |l^ 
or |l^ and ||. We will use the expression in |l^ which gives 

2m \ X 



du f ^ 2ux^ '^^^ 
x^ - 




zA\h^'--^]^ (84) 



where Cn = {u{x, t), 0+, u{x, t)). The endpoints u{x, t) of the integration must be chosen such that 
t^~" times the integral gives a solution Zi,^^{x) of the Bessel equation. This requires thatj^ 



x^ u 



2x d 



du 



x" 




= (85) 



at the endpoints of the integration contour. This condition can be satisfied for the Bessel functions 
Zi, — Jvi Id for endpoints u{x,t) = xt/2 in Cr = provided Re;/ > —1. The condition cannot be 
satisfied for other choices of the Bessel function Z^. 

We can easily show that the Riemann operator defined by ^ satisfies the product rule P^P^ = 
p^+M provided we choose the endpoints in the integrations properly. It is convenient in this to 
assume that Re A < and Re ;/ < 0, conditions which can be attained using The contour 
integrals can then be converted into ordinary integrals. The action of the group operator e^"^- on 
the integrand in ^ changes x^ to x^ — 2vx/t^ but does not affect x/t since P-{x/tY = 0. As a 
result, the parameters u and v appear only in the sum u + w as required by the operator relation 
p-iiP- p-^P- — g-(u+t))P_ _ then straightforward to show that the double integral can be reduced 
to the product of a beta function and an integral of the form in |8j, and gives a solution of the Bessel 
equation equal to f^^^^f^ Z^^x-n{z), provided the endpoints in the successive integrations are taken 
as vo{x,t,u) = xt/2 — u, uo{x,t) = xt/2. The sublety is that the endpoint of the first integration 
depends on the variable in the second. The result gives an example of the formal relation in |l^ 
which generalizes the product rule for Riemann fractional integrals. 

A change of the integration variable to v = 2u/xt converts ^ to the simpler form 



Alternatively, for Re < 0, we can collapse the integration contour in ^ and change to variables 
z = x'^, = x^ ~ 2ux/t to put the result in the form of a standard Riemann fractional integral, 

= (87) 

Re/x < 0, Rez/ > —1. This reproduces j^, 13.1(63) and 13.1(83) for Z,^ ~ J^, and Zy = when 
is replaced by — /i in accord with the convention used there. 

For Zi, = Yy, Ky, the action of the Bessel operator in the variable x = ^/z on the function 



defined by |86| leaves a term proportional to a; ^. The result is an inhomogeneous Bessel 



^The integrand vanishes for u = xt/2, suggesting that value for u{x,i). With thaL,assumed, the precise condition 
for a solution of Bessel's equation follows by scaling the integration variable as in to eliminate x and t from the 
limits of integration, applying the relevant operator, and then undoing the scaling in the resulting condition. 
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equation with a solution which involves a sum of a function J^-^ or /^-^ and the Lommcl functions 
s_i/_^+i^iy_p(a;) ||] §10.7. The fractional integral 13.1(73) is of this type. 

An analysis similar to that above shows that the integral in taken on a Riemann-type contour 
with endpoints at w = —x/2t satisfies an inhomogeneous rather than homogeneous Bessel equation 
of order v -\- The general solution involves Lommel functions, and there is no Riemann definition 
for acting on Bessel functions alone. 

3.4.2 Riemann-type integral representations for Bessel functions 

The operator relation P^t^ Z^{x) = t'^^'^Z^^^{x) immediately gives integral representations for 
and We will start with the functions of order >^ — 

2 cosx 2 coshx 

J_i/2(^) = ^^, I-iMx)^^^. (88) 

Choosing ^ = — A — ^ then gives 

Mx) = 1 -(^+^)r(-A + l)^ 



[V] / dvv^^^(\ - vY^ cos(x\/l - v) (89) 



for general A, or, replacing u by 1 — i 



2 



J^{x) = ' I dt{l-t^Y~i cosxt (90) 

«7rr(2) J(0,l+,0) 



V^r(A + 



Tt(J) / dt{l-t^Y-i cosxt (91) 
o) ^2/ Jo 



for ReA > -i. The first is a standard Poisson-typc integral representation for Jx{x) 3.3(2). 
The second gives the generalization 6.1(6). 

— A— i 1 

Similarly, from P_ '^t^2/_i (.t) = t^I\{x), 



hi.) - i^-<^^^^r(-A + i) 



^ i^Y f dvv^-^l - u)"^ cosh(a;Vl - v) (92) 
T(- — X) f 

/ dt {l-t^^-i cosh xt (93) 
»7rr(2) J(o,i+,o) 

TT o / dt (1 - t'^Y^^ cosh xt, ReA>--. (94) 
2) ^2/ Jo 2 



V^r(A 



4 Summary 

Many of the properties of the special functions arise from their connection to Lie groups 0, ||, |^. 
Their differential recurrence relations, for example, reflect the action of particular multivariable 
operators D in the associated Lie algebra, the so-called stepping operators, on the the functions in 
the relevant class. We have given general definitions of fractional operators in the context of 
Lie theory, and explored their formal properties. Our Weyl- and Riemann-type fractional operators 
generalize the single- variable Weyl and Riemann fractional integrals W-\ and 0], Chap. 13. 
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The operators 13^ change the indices on the special functions by fractional displacements related to 
A, and provide useful connections between functions in different realizations of the Lie algebra. 

We have illustrated the usefulness of the fractional operators in the case of the Euclidean group 
E(2) and the Bessel functions, and find that they contribute to a coherent overall picture of many 
relations among the Bessel functions as interpreted in the group context. For example, the formal re- 
lations P±.t^Zfj_{x) = t^^^Zf^±\{x) give the integral relations connecting Bessel functions of different 
orders. When reduced to the single variable x, these generalize known fractional integral relations. 
Used with simple choices of fi and A, with an elementary function, they lead immediately to the 
standard integral representations for the various Bessel functions, representations which are derived 
in 1^ from quite different starting points using different methods. In addition, the action of the 
elements e~"^± on the functions gives generating functions for the Bessel functions (the 

Lommel expansions), while the Bessel equation itself is the statement that the Casimir operator 
P+P- and the rotation operator iJ^ have fixed values —1 and /i. 

The applications of the fractional group operators will be extended in following paper to the 
associated Legendre functions in the somewhat more complicated case of S0(2,l) and its conformal 
extension. 
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